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Abstract 

A strongly-interacting fermion chain with supersymmetry on the lattice and open 
boundary conditions is analysed. The local coupling constants of the model are 
staggered, and the properties of the ground states as a function of the staggering 
parameter are examined. In particular, a connection between certain ground-state 
components and solutions of non-linear recursion relations associated with the Painleve 
VI equation is conjectured. Moreover, various local occupation probabilities in the 
ground state have the so-called scale-free property, and allow for an exact resummation 
in the limit of infinite system size. 

1 Introduction 

The understanding of strongly-interacting many-particle systems is one of the prime chal- 
lenges of statistical mechanics and condensed matter physics. As realistic models are a 
posteriori notoriously difficult to handle analytically, one resorts often to simplified and 
mathematically more accessible descriptions which preserve essential physical features. The 
analytical accessibility is usually due to some underlying symmetry. In this work we study 
supersymmetry, applied to models of strongly-repulsive itinerant fcrmions. 

More precisely, we consider the class of M.^ models which were introduced in [T] [2J. 
They describe spinless fermions on arbitrary lattices with the exclusion rule that connected 
particle clusters contain at most I particles. The models have a built-in supersymmetry: 
the corresponding generators - the supercharges - add or remove single particles from the 
system while preserving the exclusion constraints. We focus on £ = 1 which corresponds to 
fermions with nearest-neighbour exclusion. To date, this is certainly the best studied case, 
and presents all features of a realistic physical model. The analysis of the one-dimensional 
Ali-chain with periodic and open boundary conditions revealed interesting connections 
with enumerative combinatorics: in suitable normalisation, the ground state components are 
given by integers which enumerate alternating sign matrix number with various symmetries 
[21 [3], as are known from the properties of the XXZ spin chain at A = —1/2 (see for 
example [3J [5J |5J [7] ) . Indeed, a mapping between the models can be established, what led 
to discovery of a hidden supersymmetry in the spin chain [5]. In two dimensions, the model 
presents remarkable features [TD] such as the so-called 'superfrustration', a phase with 
extensive ground-state entropy [TT] , and direct relations to rhombus tilings [121 K31 E] ■ For 
a comprehensive introduction we refer to |15j . 
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The continuum limit of the M.\ fermion chain is described by a superconformal field 
theory with central charge c = 1 (a free boson with a special compactification radius) |16j . In 
this sense, it is quantum critical. It was understood in [T7J [TH] that an off-critical extension 
can be achieved through the introduction of a staggering parameter on the lattice. The 
staggering is 3-periodic, and closely related to the so-called 3-rule appearing in cohomology 
considerations about the model's ground states, and therefore these works focused on 
periodic lattices with the number of sites being a multiple of three. The ground states are 
polynomials in the staggering parameter. In particular, it was pointed out that several of 
its polynomial components appear also in the exact ground states of the XYZ spin chain 
along a special line of couplings, and are solutions to a non-linear recursion relations related 
to the Painleve VI differential equation jTHl HH1 HI] ■ Further evidence for the connection 
between the staggered A4i model and the special XYZ chain was given [223 : m fact, their 
Hamiltonians have coinciding spectra in several momentum subsectors. 

In this paper, we consider the staggered Mi model on a one-dimensional lattice with 
open boundary conditions. This problem has already been addressed in [23] , where several 
properties of the ground states were investigated through perturbation theory around 
trivially solvable points. In particular, the presence of kinks separating two types of ordered 
states was pointed out. The main purpose of the present article to expand upon these 
results, and show that the finite-size ground states have a number of remarkable non- 
perturbative properties. In fact, like in the case of periodic systems we conjecture a relation 
to polynomials pertaining to hierarchies of integrable equations, this time however in a 
more refined form. In particular, we point out a variety of sum rules for the square norm 
of the ground state vectors. Moreover, we exploit the so-called scale-free property of local 
observables discovered in [T7J [T5]: it states that their Taylor/perturbation expansion in the 
staggering parameter around trivial points can be summed, and thus yield a conjectured 
exact and non-perturbative form for infinite systems. 

This article is organised as follows. In section [2] we introduce the details of the Ai\ 
model, and describe different interesting variants for the staggering and number of sites. 
These variants are investigated in sections [3] and [4] we conjecture a relation between certain 
ground-state components of the models and solutions of the Painleve VI-Hirota equations, 
and provide evidence for the scale-free property of various quantities. We present our 
conclusions in section [5] 

2 The model and methods 

In this section we recall the general definition of the staggered M 1 model for fermions with 
hard-core exclusion. Moreover, we explain an iterative approach to the determination of 
its ground states, exploiting the fact that it is a polynomial in the staggering parameter. 
The key polynomials appear to be deeply related to the theory of the Painleve VI equation, 
and were first discovered in the context of a particular XYZ chain []1|] EH] . We recall briefly 
their main features. 

2.1 Definition of the staggered M.\ model 

We consider a one-dimensional chain with N lattice sites which we label by integers 
j = 1,2, . . . , N . The model describes spinless fermions living on this chain. These are 
created and annihilated by operators Cj , cj which obey the usual anti-commutation rules 

{c,-, Cj} — {cj, c^} = and {c^, ct} = The particle configurations are constrained by a 
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nearest-neighbour exclusion: at most one of two adjacent sites may be occupied. We will 
represent an occupied site by 1, and an empty site by 0. For example, a configuration for a 
six-site chain with particles on sites 2 and 6 is given by a = 010001. The hard-core exclusion 
forbids therefore pairs like • • • 11 • • • . We will sometimes abbreviate the empty configuration 
by = 00 • • • 0. We shall work with free boundary conditions what is equivalent to add two 
inaccessible but empty sites at j = and j = N + 1. 

Our aim is to study a Hamiltonian generated by the two supercharges 

N N 

q = E a a £ t = E A i4 

where dj — Pj—iCjPj+i, with the projector Pj = 1 — CjCj, are dressed fermion operators 
which respect the hard-core exclusion. As the sites and N + 1 are always empty, we 
set Pq = Pn+i = 1- The X 3 ■, j — 1, . . . , N are arbitrary complex numbers which we call 
sometimes coupling constants. The supercharges are nilpotent for any choice of the X/s: 
Q 2 = {Q ] f = 0. The Hamiltonian is given by H = {Q, Q^}. In terms of the fermions the 
H it can be written as 

N-l N 

H commutes with the fermion number operator F — Y2j=i c ] c j> anc ^ °f course the super- 
charges. The algebraic relations of the quadruple H, Q,Q\F yields the well-known M = 2 
supersymmetry algebra |24j . 

In this article we concentrate on the case where the number of sites is of the form 
N = 3n or N — 3n — 1. For both choices and non-zero coupling constants the model 
has a single zero-energy ground state in the subspace of n particles as can be shown by 
cohomology arguments [14] , Conversely, for N = 3n — 2 no such zero-energy ground 
state exists: supersymmetry is spontaneously broken. The proof of this statement uses a 
one-to-one correspondence between linearly independent ground states and the elements of 
the quotient space (or cohomology) F)q — ker Q/im Q. Its dimension coincides with the 
dimension of the ground state space. In order to determine F)q, one proceeds into two steps. 
The lattice is divided into two disjoint parts, call them Si and S2. As the supercharges 
are linear superpositions of fermion annihilators/creators, we divide Q — Q\ + Q2 so that 
Q1/2 a cts only on S\j2- First, one evaluates $Jq 1 - Second one acts with Q2 within this 
space in order to find f)Q, applying the so-called tic-tac-toe lemma. Even though the details 
of this procedure are not important for our considerations, we wish to point out that the 
division into sublattices uses a fundamental feature of the M. \ model, namely an underlying 
periodicity with period 3. Indeed, in the existence proof Si = {2,5,8,...} is chosen to 
contain every third site, and S2 all the remaining sites. This inherent "3-rule" led to the idea 
to choose the coupling constants Xj of the supercharges from a 3-periodic pattern [T71 [TB] , 
namely Ai = y, A2 = 1, A3 = y, . . . , and investigate their properties as a function of the real 
parameter y. For periodic chains with N — 3n sites, it does not matter which one of the 
three coupling constants is scaled to one. For open chains however, various choices lead to 
different situations. We choose to investigate two one-parameter staggerings: 

(I) X 3p - 2 = y, X 3p -i = 1, A 3p = y, (1) 
(II) A 3p _2 = y, A 3p _i = y, X 3p = 1, (2) 
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where p = 1, . . . , n. Notice that the choice (I) is invariant under a parity operation P which 
reverses the order of all sites^ For this choice, the Hamiltonian commutes with the parity 
operator [H, P] = 0, and eigenstates may be classified as parity-even and parity-odd. Under 
parity the staggering pattern (y, y, 1) of (II) is mapped to (II') (1, y, y), whose properties 
can readily be deduced from (II), and will therefore not be considered separately. As y = 1 
all the cases reduce to the quantum critical model studied in [3J. Here, we focus on general 
values for y. 

2.2 Polynomiality and diagonalisation methods 

Our aim is to find the properties of the zero-energy ground states. It is natural to choose to 
do this in the occupation number basis. Therefore we expand the single solution to solving 
H(y)\*f?(y)) =0 according to 

l*(y)> = Z>«(v)|a>. 

a 

Here a runs over all admissible configurations on N sites with n particles, where N = 3n 
or N = 3n — 1 according to the case under consideration. For both choices (I) and (II) of 
the staggering, the supercharges are linear in the coupling constant Q = Q + yQi, and 
therefore the Hamiltonian is a quadratic polynomial in y: H(y) = Hq + yH± + y 2 H 2 . This 
implies that we can choose the normalisation for \^(y)} in such a way that all components 
ipa{y) are polynomials. As this choice might still contain some redundancies, we impose 
that for any given configuration a there is at least one a' such that the polynomials ip a (y) 
and tp a '(y) have no common factors (i.e they are coprime). This fixes the ground states up 
to an overall c- number normalisation which will be adjusted accordingly. 

The actual construction of |^(y)) is based on the following observation. Since all ip a {y) 
are polynomials, we can expand 

i*G/)>=xyi*<>' 

for a certain degree d, depending on the lattice size N. In all the cases studied below, we 
provide the explicit form of the zero-order term |^o) which is obtained in the limit y — >• 0. 
Of course, this is a well defined state, although the zero energy space can be degenerate at 
y = 0. The above decomposition can also be written as 

|*(l/))=m(tf)|* > + !*(»)>, 

where m(y) = 1 + Yli=i m t f° r a certain degree d' , and (^ol 1 ! 7 ) = 0. In all cases, we 
provide a conjecture for the polynomial m(y) (again, depending on the lattice size). This 
means that the identification of the non-degenerate state obeying H \ ^(y)) = 0, for generic 
y, will reduce to the solution of 

H + Hi + H 2 |* £ _ 2 ) = 0, 

where the ambiguity in the component of ker Hq is solved in terms of m(y). Indeed, by 
consistency we have 

H 1 + H 2 \V t _ 2 ) _L ker H , 

J We leave aside the question of possible signs through the reordering of fermions. 
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and we can write 



This simple algorithm can be applied starting with = | \I/ 2 ^ = and, by consistency, 

terminates as soon as \§ 1) — l^-i) = for a certain I — d. Here H^ 1 is the inverse on 
the orthogonal complement of keriJo- Its computation is trivial, since Hq is diagonal in 
the occupation number basis. Moreover, we exploit sparse-matrix linear algebra in order to 
apply Hi and Hi. 

2.3 Relation to the XYZ chain and Hirota equations 

It was pointed out in jTTl [TB] that the ground states of the staggered fermion chain with 
periodic boundary conditions and N = 3n sites present some striking similarities to the 
ground states of the XYZ Hamiltonian 

I L 

k — 1 a—x,y,z 

Here a a , a = x,y, z denote the usual Pauli matrices. We choose periodic boundary conditions, 
and restrict the coupling constants J a to the the line 

For L = 2n + 1 sites this chain has exactly two translationally-invariant ground states |$±) 
at the energy E n = —L((^ 2 + 3)/4 [25]. These ground states are related by the reversal 
of all spins, and therefore it is sufficient to study one of them, say |$_) which contains 
the configuration || ■ • • J,) but has zero overlap with |f ■ ■ ■ t)- Bazhanov and Mangazeev 
showed that under suitable normalisation the component in front of the completely polarised 
configuration 44- " " 4- can b e written as 

^u-i(C) = C in+1) s n (C 2 ) 

where the s n (z), n S C are polynomials in z of degree \n 2 /4\ which solve the non-linear 
differential recursion relation 

2z{z - l)(9z - l) 2 (lns„0))" + 2(3z - l) 2 (9z - l)(lns„(z))' 

+ 8(2n + lf ^%^ - (4(3n + l)(3n + 2) + (9z - l)n(5n + 3)) 

= 0, (3) 

with initial conditions so(w) = si(w) = 1. This recursion is a special case of the Hirota 
equations satisfied by a tau-function hierarchy associated to the Painleve VI differential 
equation [251 H5]. 

These polynomials appear also in the two zero-energy ground states of the staggered 
fermion chain at N = 3n sites with periodic boundary conditions. The precise relation is as 
follows. The ground states can be chosen to be parity eigenstates. Write them as |^' ± (?;)) 
where the label ± corresponds to the parity eigenvalue. Upon appropriate normalisation we 
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find that the components of the configuration 010010 • • ■ 010 are related to the s-polynomials 
according to 

7 + , s _ C L( - 1)2/2J «n-l(r 2 ) 

VOIO-OIOW - 2 („-l)(n-2)/2 s^^O) ' 

I- f , - ^}^_S-r^) 

VOIO-OIOW - 2 („-l)(„+2)/2 s _ n (Q) ' 

where the variables y and £ are related through the quadratic equation 
C 2 = l + 8y 2 . 

This hints strongly at a profound relation between the two models, but a precise mapping 
is so far only known in the critical case y — 1, £ = 3. 

2.4 Factorisation of Hirota polynomials 

In this work, we will point out a connection of the open staggered chains to the Hirota 
polynomials. It is more subtle than for periodic boundary conditions, and related to 
factorisation properties of the s n (z). These were conjectured in [50], and here we quote this 
conjecture: 

Conjecture 2.1 For all n G Z the ratios s n (z)/s n (0) are polynomials with integer coef- 
ficients. Also, for all n > — 1, s„(0) = 1. For any n they factorise into two polynomials 
according to the following scheme: 

1. For n = 2k + 1, k G Z, we have 

S2/C+1O 2 ) = s 2 k+i(0)p k (w)pk(-w), pfe(0) = 1. (4) 

where the Pk{w) are polynomials with degpk(w) = k{k + 1) and integer coefficients. 
We have the transformation property 

1 + 6W \ I 1 — W 



2. For n = 2fc, fceZ, i/ie polynomial factorises according to 

S2k(w 2 )^c k (l + 3w) k(k+ ^p^k-i(^ I ^-)qk-i(w), g fe (0) = l. (5) 

\1 + 3w/ 

Cfe = 2 _ ' £ (' c+2 ) if k > 0, and 2~ fc2 (2/3) 2fe+1 otherwise. The polynomials Qk(w) 
have deggfe(w) = fc(fc + 1) and positive integer coefficients. Moreover, they are even 
and obey the transformation rule 

, , /l + 3w\ fe(fc+1) / 

A list of the first polynomials is given in the appendix. In the next sections, we will 
present several conjectures about the ground states of the open staggered fermion chain. 
In particular, we claim that the polynomials Pk(w) and Qk(w) appear in various disguises, 
not only as components in the ground state vectors but also in their square norms, and 
projections on the ground states at weak coupling y = and strong coupling y — >• 00. 
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3 Ground-state properties of the staggered chain with 
N = 3n sites 



We start with staggered open chain for the two choices (I) and (II) for N = 3n sites. 
In both cases, we first determine systematically the ground states at y = and y — s- oo 
(weak and strong coupling limits). Next, we present results at finite coupling y from exact 
diagonalisation of small systems which we use as a basis for our conjectures. Finally, we 
study site occupation probabilities near the boundaries and in the bulk, point out scale-free 
expansions in around the asymptotic points, and use them to conjecture exact expressions 
in the limit of infinite systems. 



3.1 Staggering I 

3.1.1 Asymptotic ground states 

We begin with the derivation of the y = ground state, which turns out to be a kink state 
as was shown in \2'6\ by means of perturbation theory. Here we prefer rather to exploit the 
properties of the supercharges Q, alone. The considerations are very reminiscent of the 
3-rule, and use polynomiality/analyticity in y. Moreover, they will prove to work also for 



staggering II (see section 4.2) 



At y = 0, we have the supercharge Qq = Y^j=i daj-i- As it acts on isolated sites, it is 
clear that any state annihilated by this operator must not contain any particles on the sites 
3j — 1 with j = 1, ... ,n. For Qq to annihilate such a state however, there must a least be 
one particle adjacent to such a site. This requirement is fullfilled by all states of the form 

|m) = 1100100- •• 100001001 •••001), m = 0,...,n 

3m 3(n— m) 

These are kink states in the sense that they interpolate between two ordered configurations 
100100 • • • and • • • 001001. All of them are admissible ground states at y = 0, and therefore 
E = has an (n + l)-fold degeneracy. The general solution is a superposition \^(y = 0)) = 
\^k) = Sm=o °m|w)- In order to determine the relevant combination we demand continuity 
with y ^ 0, where we know that only a single ground state may exist. The polynomial 
character of the ground-state wave function allows to write \^(y)) — \ 1 &K)+y\ 1 &i) + - ■ ■ Acting 
with the supercharges leads to the equations Qo\^i)+Qi\^k) = and Ql\^ i) + Q\\^ k) = 
at first order in y. We claim that the second equation is already sufficient to determine 
\^k), even without knowing the structure of To show this, we examine the action of 

Q\ on the states |m). Let us introduce the defect states 

\<f>m) = | 100- - - 100 101 001 ■ • ■001 ), m = 0,...,n-l 

3m 3(71— m— 1) 



with Ti+l particles. Then we have 
Q\\m) = (~iy 




m = 0, 
m = n, 
otherwise. 



None of the states \(f> m ) contains any particles on the sites 3j — 1, and therefore none of them 
can be obtained by acting with Qq on any other state having n particles. So we must have 
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QqI^i) = 0. This implies the additional requirement that QII^k) = 0. More explicitly, we 
find 

n-l 

^(-l) m (a m -a m+1 )|0 m ) =0. 

As the \4>m) are linearly independent, we conclude that ao = <2i = ai = ■ ■ ■ = a n . This 
determines already the ground state up to an overall normalisation because for y > it is 
unique. Setting the normalisation to unity we thus find the ground state 

n 
m=0 

which is a uniform superposition of kinks. This is the result of |23j . 

As y — > oo, the ground state has to be annihilated by Qi and its adjoint Q\. Let us write 
this in an suggestive way: Q\ — di + Y^j—i (^3j + ^3j'+i) + <^3n- From this it is clear that 
any state annihilated by both these supercharges must not have a particle on the boundary 
sites 1 and 3n, but a particle on their neighbours 2 and 3n — 1 respectively. Because of the 
hard-core exclusion, the situation obtained in this way is equivalent to the initial one with 
six less sites (unless one starts with only three sites where the problem is trivial) , and one 
can proceed by iteration. The ground state for y — > oo will therefore be proportional to the 
polarised state 010010 • • • 010: 

n 

hM = E4j-iio>. (6) 

3.1.2 Finite coupling 

Now we turn to finite values for the coupling constant y. We list some properties, based on 
observations from exact diagonalisation of the Hamiltonian for small system up to N = 24 
sites (n = 8 particles). In order to give a simple illustration, we display here the ground 
state components for six sites. As the Hamiltonian commutes with parity it is sufficient to 
indicate the representatives 

101000 -y 2 100100 l + 3y 2 

010100 -y{\+Ay 2 ) 100010 -2y (l + 2y 2 ) 

010010 y 2 (3 + 8y 2 ) 100001 1 + 2y 2 

The ground state wave function is in fact even under the parity operation. It is not difficult 
to see that this holds for any n: the asymptotic ground states at y — > and y — > oo are 
parity-even. Since the eigenvalues of P are discrete this must also hold at finite coupling. 
From our example we observe that all components have definite transformation behaviour 
under y — > — y: if S a denotes the total number of particles on the sites 2,5,8,. . . , 3n — 1 
in configuration a then we have ip a (—y) = (—l-) Sa 4>a(y)- This is related to the fact that 
the unitary operator U = expi7r^" =1 c^^^sj-i has the property UH(y)U^ 1 = —H{—y). 
Taking into account the uniqueness of the ground state and the boundary condition at 
y = 0, we conclude \^{-y)) = U\^{y)). 

Let us now have a closer look at certain components. To motivative this, we compute for 
our example of six sites the square norm of the ground state vector, as well as its projections 
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on the asymptotic ground states at y = and y — > oo: 

||*(y)|| 2 = (3 + 8y 2 ) (l + 6y 2 + lly 4 + 8y 6 ), 
(* P |*(y)) = y 2 (3 + 8y 2 ), (*at|*(»)) = 3 + 8y 2 . 

The polynomial 3 + 8y 2 appears thus in multiple locations. It is one of the two factors of 
the square norm. The other factor can in fact be written in terms of the polynomial q2{w) 
related to the Hirota equation, after a suitable change of variables to £ = yT + 8y 2 . The 
systematic examination of the ground states up to n = 8 leads to the following 

Conjecture 3.1 All ground state components are polynomials with integer coefficients and 
definite parity under y — » —y. The projections on the kink state and the polarised state 
differ only by a monomial 

(* P \H>(y)) = y n (* K my)) 

Here (9k\^{u)) zs an even polynomial in y of degree n(n — 1) with integer coefficients. The 
square norm factorises according to 

/ f \ n(n-l) 

\My)\\ 2 = (*KMy))x [|J ^(C 1 ), C 2 = i + 8y 2 , 

where q n (w) is the polynomial implicitly defined in 
The general form of the ground state is 

n 

|¥(l/)> = £>k(lO l*> + l*>. (fc|*) = 0, 

fc=0 

with certain polynomials mfc(y), and \k) the individual kink states defined in section. We 
normalise them according to rn^iy) = 1 + 0{y 2 ). They are fixed by the requirement that 
the iterative iterative algorithm for the construction of \^{y)) terminates. 

We have not found a relation between the projections on the kink states, and the 
polynomials coming from the Hirota equation ([3]). Here we list t n (y) = (^k\^(u)) f° r 
n = 1, . . . , 5: 

h(y) -2, 

t 2 (y) = 3 + 82/ 2 , 

i 3 (2/)=2(l+4y 2 ) (2 + 7y 2 + 8y 4 ), 

f 4 (y) = 5 + 72y 2 + 420y 4 + 1300y 6 + 2304y 8 + 2304y 10 + 1024y 12 , 
t 5 (y) = 2(3 + 70y 2 + 720y 4 + 4332y 6 + 17020y 8 + 46053y 10 + 87840y 12 
+ 118080y 14 + 108288y 16 + 61440y 18 + 16384y 20 ), 

3.1.3 Densities 

We probe the system by examination of the probability to find a particle at site j: 
p«(j/) = <*(y)| C tc i |*(y))/||*(y)|| 2 . 
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In particular, at the boundary j = 1 we know the result for any n at the critical point y = 1 

01 

,!»(y = D= " (10n + 11) 



2(4n + 3)(4n + 5) 



Therefore, in the limit of large systems n — > oo we find the density 5/16. 

For arbitrary y the densities are given by ratios of polynomials at finite n. While we do 
not have closed-form expressions for these we examine their series expansions around the 
points y — and y — oo, by following the strategy of |18j . There it was pointed out that 
these series expansions have a so-called scale-free property: the first ~ N coefficients do not 
depend on N. Let us examine this for the boundary density at site j = 1 near the strong 
coupling point y — oo. Up to n — 6 the Taylor expansions of Ay 2 Pn\y) are given by 

4y 2 Pl (y) = 1-1+^-^+^ + ... 

An 2 n (n\ = 1 4- -J- - 39 | 321 _ 1735 , 6209 5671 

P2\y) i T 8y 2 64 , y 4 -r 512y 6 4096j/ 8 32768y 10 262144; 



•y P2\y) L "r s„2 64y 4 ~r 512y 6 4096y 8 32768y 10 262144t, 12 

A y 2 ( \ _ 1 i _i_ i _3 417_ , 5347 39505 , 201147 

l J 3\y) 1 i 8y 2 i 64 „4 512 e i 4096i/ 8 32768y 10 ~ 262144t/ 12 

4,,2„ ,C 7 A _ 1 i JL i _i 3 _ 2403 , 42355 , 101479 

^i/ ' 8y 2 64i/ 4 T 128y 6 2048y 8 ^ 16384y 10 ^ 32768?, 12 

4v 2 n^(v) = 1-1 ± U -I ? , 7 29249 , 659389 

^y l J 5\y) ~ 8j/ 2 ~ 64j/ 4 ' 128y 6 ~ b!2y s 16384y 10 ^ 131072)/ 

Al .2 ( \ I I 1 I a _i 3 _|_ __7 9 185329 

^y P6\y) ± T 8y 2 T 64y 4 T I28y 6 > 512y 8 1024i/ 10 65536J/ 12 



The bold-printed terms correspond to the stable part of the series, i.e. they are not changed 
as n increases. The stable part of the series matches the series expansions 



l + 12y 2 -8y 4 + 8y (y 2 - if 2 

— i y ^ ±- 



We conjecture that this is indeed the thermodynamic limit of p n (y), however the singularity 
at y = 1 with critical exponent 3/2 indicates that this can only hold for y > 1. This conjecture 
is in fact supported by the observation that this function reproduces the correct asymptotic 
value Poo(y — 1) = 5/16 at the critical point. 

Around y = we observe that the series expansions of boundary density does not have 
the scale-free property. We attribute this to the non-local nature of the kink state, as 
opposed to the very local and ordered character of the polarised state at y — >• oo. The first 
coefficients match the pattern 

fiw \ 11 ( \ 71 + 3 2 

Pn(v) = —rr \ 1 -—rTV +■ 



Yet we are able to guess the limit as n — > oo. To this end let us anticipate a little. In the 
next section we study the chain with staggering (I) and the number of sites N = in — 1 . As 
it will turn out, that case has a simple ordered state at y = 0, and allows to evaluate the 
corresponding boundary density pn\y) near y = from a scale-free expansion. In the limit 
of infinite system size n — > oo the expressions in these two cases should match as the other 
boundary is very far away. Thus we conjecture that also in the present case N = 3n the 
limit n — > oo yields 

Poc w 32y2 , u - 
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Next, we discuss the centre of the chain. In the limit of large systems, we expect the 
average values of local densities to measured in the bulk to be insensible to the boundary. 
Therefore it is reasonable to speculate that some of them will be related to local densities of 
the closed chain with periodic boundary conditions in the thermodynamic limit. We probe 
the density at the site j c — 3n/2 — 1 for even n, and j c = (3n + l)/2 for odd n. Despite 
the presence of the kink state in the weak-coupling regime we find scale- free expansions. 
Extracting the pattern for the stable part, we conjecture the following expressions 

( i+2y*-y/l+8y* , 

I i + 8 y 2 , y^x 

These correspond to the staggered densities of the closed chain in the thermodynamic limit, 
for the parity-even ground state |18j . The parity sector is expected because in the present 
case the open chain ground state is parity-even, too, as shown above. 



3.2 Staggering II 

3.2.1 Asymptotic ground states 

Let us first evaluate the ground states in the limits y — and y — > oo. At y = the 
ground state has a simple structure as opposed to the kink state of case (I). We have 
Q(y = 0) = Qo = 52"=i d>3j. For Qq to annihilate a state, there must not be any particle on 

the sites 3j, j = 1, . . . , n. Conversely, if Qq is to annihilate such a state, each site 3j must 
have at least one neighbour occupied by a particle. All sites 3j have two neighbours except 
for the last one 3n. Hence there must be a particle on site 3n — 1, and by iteration it is 
not difficult to see that this implies that all sites 3j — 1 have to be occupied. Hence a state 
annihilated by both Qo and Q is given by 010010010 ■ ■ ■ : 

n 

|* P ) = |*(y = 0)> = ni--i|0>. 

i=i 

Second, as y — > oo the ground state must be annihilated by Q\ = Ylj=i ^3j-2 + ^3i-i> 
and its Hermitian conjugate Q\. In the sector with N — 3n sites and n particles, Q[ 
annihilates any configuration that contains a particle in each on the n pairs of sites 
(3j — 2,3j — 1), j = 1, ... ,n. As these are disjoint, we can solve the problem for Q\ 
separately for each such pair. Consider thus the sites 1 and 2, with adjacent sites and 3 
being empty. The only way to produce a state that is annihilated by d\ + d,2 with a single 
particle is given by the singlet ••• 10 — 01 ••• . We conclude that the ground state 
becomes a tensor product of such singlets 

n 

i*s> = II(4 J -- a -4i-i)|o>. 

3=1 



3.2.2 Finite coupling 

Now we consider general values for y. We shall normalise the polynomial coefficients in the 
expansion \^{y)) = *Pa{y)\a) such that the zero-order term for a — 010010010 • • • is +1. 
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As an example, we present the ground state for N — 6 sites in the sector with n = 2 
particles. We find the amplitudes 

101000 -y 5 100100 2y 4 (l + y 2 ) 

010100 - y 2 (1 + 2y 2 + 2y 4 ) 100010 - y 2 (l + y 2 ) (l + 2y 2 ) 

010010 1 + Ay 2 + 4y 4 + 2y 6 001010 - y (l + 3y 2 + y 4 ) 

100001 y 3 (1 + 2y 2 ) 010001 - y (l + y 2 ) (l + 2y 2 ) 

001001 y 2 (l + 3y 2 ) 000101 - y 3 

We see that all components have integer coefficients and definite parity under y — > —y. The 
second property can be justified in a way which is similar to the one for staggering (I). The 
projections on the asymptotic ground states are given by 

(*p|*(»)) = l + 4y 2 + 4y 4 + 2y 6 , 

(*s|*(l/)> = l + 6y 2 + lly 4 + 8y 6 . 

We find these two polynomials as factors in the square norm. Indeed, for six sites 

||*(y)|| 2 = (l + 4y 2 + 4y 4 + 2y 6 ) (l + 6y 2 + lly 4 + 8y 6 ) . 

We have verified these properties up to n — 8 particles (N = 24 sites), and identified the 
projections on the asymptotic ground states in terms of polynomials associated with the 
Hirota equation after suitable variable transformations. Thus, we formulate the following 

Conjecture 3.2 With the above normalisations all ground-state components are polyno- 
mials with integer coefficients and definite parity under y — > —y. For all n we have the 
projection 

(*p|*(»)> - {^) nin+1 \ n (*^f) , C 2 = 1 + Sy 2 . 

As I'J/p) is a single- component state, we can generate ground states for general n with the 
help of the iterative algorithm and m(y) = p\*f> (y)) . Likewise, the projection on the 
singlet state is given by 

/A n(n+l) 

<*s|*(y)>= M q^C 1 ), C 2 = l + 8y 2 . 



2 / 

The square norm factorises into the product of these two projections 

||v]/(y)|| 2 H*p|*(y))(*sl*(y)}- 

3.2.3 Densities 

Let us again examine the probability to find a particle at the first site. For n — 1 , . . . , 4 we 
find that it is scale-free: 

y~*Pi\y) = 1 - 3y 2 + 7y 4 - 15y 6 + 31y 8 - 63y 10 + 127y 12 + . . . 
y^P2 ] {v) = 1 - 3y 2 + 12y 4 - 56y 6 + 246y 8 - 1002y 10 + 3884y 12 + . . . 
y~ 4 Ps\y) = 1 - 3y 2 + 12y 4 - 56y 6 + 288y 8 - 1584y 10 + 8723y 12 + . . . 
y~ A Pi\y) = 1 - 3y 2 + 12y 4 - 56y 6 + 288y 8 - 1584y 10 + 9152y 12 + . . . 
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We checked the scale-free property up to n — 8 particles. The stable coefficients are correctly 
reproduced by the generating function 

which gives the correct value Po^{y = 1) = 5/16 at the critical point. Around y — oo the 
quantity is scale free, too. We find 



P { i\y) 



1 3 7 15 31 

2 ~ V + Sy 1 ~ 16y6 + 32y* 
1 1 65 1083 8601 



o^tv) = — + """" ™^ 

P2 w 2 8y 2 64y 4 512y 6 4096y 8 

P ( i\y) = 



1 1 1 27 607 



2 8y 2 32y 4 16y 6 128y 8 

(i) 1 1 1 3 12179 

Pi \y) - 2 8y2 32y4 256y6 4096^8 

Again, we checked this up to n — 7 particles, and inferred the infinite system limit: 



1- V + 8y 4 -8y (y 2 - if 2 
16y 2 



(y) = — — \a 2 — — . y^ 1 



Notice that this has the same structure as the result for the case (I) staggering, and is 
compatible with the value 5/16 at y = 1. This guess implies a discontinuity for the second 
and higher derivatives at the critical point. 

As similar pattern is found by probing if a particle is present on the last site. The data 
is consistent with p^\y) — y~ 2 Pn\y) for finite n < 8. We conjecture this to hold for 
arbitrary system sizes. The expectation values in the thermodynamic limit are then readily 
deduced from our previous conjectures. 

Let us now consider the centre of the chain. We focus on the density Pn"\y) with 
j c = 3n/2 for even n and j c = 3(n + l)/2 for odd n. It turns out that this quantity is 
scale-free both in the weak-coupling and the strong-coupling regime. From the stable part 
of the corresponding series expansions we infer the expression for infinite system size: 



— i+8p — . y ^ 1 



These coincide with the expectation values for the densities in the parity-odd ground state 
for the staggered periodic chain. 



4 Ground state properties of the staggered chain with 
N = 3n — 1 sites 

In this section we consider N = 3n — 1 sites. Despite the fact that the staggering is 
'incommensurate" we will point out that a relation to the Hirota polynomials continues to be 
present. Moreover, we shall show that the kink state appears in the weak-coupling limit of 
staggering (II) whereas the asymptotic ground states for staggering (I) have a very simple 
structure. 
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4.1 Staggering I 



4.1.1 Asymptotic ground states 

Let us focus on the chain with N = 3n — 1 sites and staggering (y, l,y). We determine the 
asymptotic ground states which are considerably simpler than in the cases studied so far. 
For y — > we seek for a state that is annihilated by Q = ^3j-i and its adjoint. So we 

know that for such a state every site 3j — 1 must be empty, but have a neighbouring site 
which is occupied. The site 3n — 1 has only one neighbour and therefore we must place a 
particle at site 3n — 2. Because of the hard-core exclusion there cannot be a particle on site 
3n — 3. However, since the site 3(n — 1) — 1 must have a particle in its neighbourhood we 
are left with the same problem as before, this time with the number of sites reduced by 
three, and the number of particles by one. We can thus proceed by iteration and see that 
the ground state is given up to a multiplicative factor by 

\^(y = 0)) = |*l> = |100100 • ■ ■ 10010). 

Moreover, we see that the kink for N = 3n sites at y = basically appeared because we 
had the freedom to put a particle on the additional site 3n. 

For y — > oo the construction of the ground state is similarly straightforward. It must be 
annihilated by Q — d\ + Y^j=i ^3j + ^3j-i and adjoint. As Q acts on the isolated site 
1 we conclude that this one has to be empty but contain a particle in its neighbourhood, 
thus there must be a particle on site 2. Because of the nearest-neighbour exclusion, there 
cannot be a particle on site 3, nor is it possible to insert one there. Therefore the state 
is automatically annihilated by and its adjoint. We see that we are now left with the 
problem of placing n — 1 particles on a chain with 3(n — 1) — 1 according to the same rules 
as before. Again, we can proceed by iteration. Therefore, the only possible ground state at 
strong coupling is given up to normalisation by 

~ |*r) = 1010010- •• 01001). 

The simplicity of these asymptotic ground states will lead to simple scale-free expansions as 
we shall see. Moreover, notice that and represent exactly the two type of order 
which the kink state of section [3 . 1 . 1 1 interpolates . 

4.1.2 Finite coupling 

For TV = 5 sites, the ground state components at finite coupling are given by 

10100 -y 2 10010 l + 2y 2 

01010 -y(l + 2y 2 ) 10001 - y (l + y 2 ) 

01001 y 2 {l + 2y 2 ) 00101 - y 3 

We see that all components are polynomials with integer coefficients which have a definite 
sign for a given configuration. The sign is flipped when going from one configuration to 
another which differs from the former by a single fermion move. The parity of the component 
polynomials under y — > — y can be obtained by counting the number of particles on sites 
3 j — 1. This is related to a simple unitary transformation implementing this reflection, 
as explained previously. The square norm is ||\I>(?/)|| 2 = (l + 2y 2 ) (l + Ay 2 + Ay A + 2y 6 ). 
We have L \*!> (y)) = (1 + 2y 2 ) and (* fl |*(y)) = y 2 (l + 2y 2 ), i.e. one of the polynomials 
occurring in the square norm can be obtained as projections on asymptotic ground states. 
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Conjecture 4.1 For general n the ground-state components are polynomials with integer 
coefficients and definite parity under y — » —y. We have the projections 



/ f\ n(n-l) 

<*i|*(i/)> = (^J Qn-iiC 1 ), C 2 -l + 8y 2 

and 



= (-l/)"<*i|*(»)>. 

^4s &oi/i \^l) and \$>r) are single- component states, we can straightforwardly apply the 
iterative algorithm with m(y) = ($! L \ty (y)) to generate general ground states with arbitrary 
n. These ground states are therefore characterised by a single polynomial. The square norm 
factorises into two polynomials with integer coefficients: 

/H A" (n+1) /I -A 
||*(»)|| 2 = <*L|*(l/)>x^-pj P™(3+i)' C2 = 1 + 8y2 ' 



4.1.3 Densities 

The boundary densities at site j = 1 have scale free series expansions both at y — and 
y — > oo. From the stable parts of their series we conjecture their forms in the limit of infinite 
systems: 



l + 44y 2 -8y 4 -(l+8;/ 2 ) 3/2 ,. < i 
nW (,,)-) 32y2 > V SI 

^oo W ^ l + 12y 2 -8y 4 +8yfa 2 -l) 3 ~ 



k rep . y^ 1 - 

As anticipated earlier, the expression for y < 1 is equivalent to the one for N = 3n sites. 



4.2 Staggering II 

We now turn to staggering II with patter (y, y, 1). As it will turn out, for N = 3n — I this 
is the choice which leads to a kink state at y = 0. The derivation of its precise form turns 
out to be slightly more complicated than the N = 3n case. 



4.2.1 Asymptotic ground states 

We start with y = 0. The supercharge is given by Q = Qq = Y^j=i^j- We focus on 
n > 1. A state annihilated by Qo and Qj must not have particles on sites 3j but at least 
one particle adjacent to them. We will call this property "adjacency condition". We see 
immediately that it can only be fulfilled if there are at least n — 1 particles, and at most n. 
The second sector contains the ground states which will be compatible with non-zero values 
of y. The sites 3j divide the system naturally into n boxes of 2 sites. Here is an example 
for n = 3: 



□□ 





□□ 





□□ 



Given n particles we need put a particle in each of them. Consider the leftmost of these boxes. 
If we put a particle on site 1 then there must be particles on all sites 3j — 2, j = 1, . . . , n in 
order to satisfy the adjacency condition. However, if the particle in the leftmost box is at 
site 2, we already satisfy the adjacency condition for site 3. For the remaining n — 1 boxes 
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we are left with the same problem as before, but with the number of sites reduced by three. 
This reminds us of the kink state problem which we encountered for the staggering (y, 1, y) 
in section [3. 1.1 1 in fact, any of the following n + 1 kink states 

\m) = | piO01() ■ ■ ■ OIOIOOICX) ■ ■ ■ 10O 10), ro = 0, . . . , n - 1, (8a) 

m n—1— m 

\n) = I01O010O- ••01O01) (8b) 

is annihilated by the supercharges. We see that for 1 < m < n — 1 the state |m) contains a 
defect of type • • • 010100 • • • . We expect to ground state of our problem at y = to be a 
linear combination: 

n 

1*0) = E °™i m )' 

m=0 

where the a m are complex numbers to be determined from continuity to arbitrary values of the 
coupling. For finite y the wave function will be a polynomial = l^o) + 2/|^i) + • • • with 

relations Qo|*i) + Qi|*o) = 0, and Qj|*i) + Ql|*o) = where Qi = E"=i( d 3^2 + d 3j -_i). 
All states in \*S> ) are already annihilated by Q[. So we concentrate on the action of Qi on 
|^o) ■ This operator takes out particles from the n boxes introduced previously. Thus, it is 
useful to introduce states \m,j), j — 1, . . . ,n which correspond to the state \m) with the 
particle from the j-th box removed. This definition is somewhat redundant because of the 
identity |m — 1, m) — \m, m) for 1 < m < n. So let us write explicitly 

n n 

where the sign comes from the fermionic character of the particles. As we wish to solve the 
equation Qol^i) + Qil^o) = 0we should ask if a state \m,j) can be created by acting with 
Qo on something (i.e. if it is Qo-exact). First, for m = this is possible if j = 2, 3, . . . , n 
because only then there will be an empty site 3(j — 1) with no particle on a neighbouring 
site. The illustration 



|0,1> = |00O10O10O---10O10) 

|0,2) = |10O00O10O---10O10) 

t 

|0, 3) = I10O10O00O10O • • • 10O10) 



points out that the state |0, 1) cannot be created through the action of Qo- A very similar 
argument holds for \n, n) which cannot be obtained through Qo acting on some other state, 
whereas the \n, j), 1 < j < n can. Let us now turn to the more general 1 < m < n — 1. 
A state \m,j) can be created through the action of Qo if j 7^ to, to + 1. Again this is 
because in these cases removing particles will create an empty site 3j with no particle on a 
neighbouring site. If however j = m or j = m + 1 then there must be a particle on one of 
the sites adjacent to 3j, and so in this case a creation through action of Qo is not possible. 

Let us now write Qi l^o) = |^o) + l^o) where \^' ) collects all states which can be created 
through the action of Q on something, and I'J'q) those which cannot (in mathematical 
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parlance, \^' Q ) is Qo-exact, while \^q) is not). By definition, we must have |^g) = 0. More 
explicitly, we find 

l*o) = = E ((-l) m_1 am|m,m) + (-l) m a m \m, m + 1)) 

771=1 

+ a o |0,l) + (-l)"- 1 a„|n,n). 
Using now the rule \m — 1, m) = |m, m) , this can be reduced to 

n 

53(-l) n, - 1 (a in _i + a m )|m,m) =0. 

m— 1 

Because all states in this sum are linearly independent, we find that a = —a.i = a 2 = 

— 03 = We fix the normalisation by setting ao = 1. This completes the construction of 

the ground state at y = and shows that it is again a superposition of kink states, in this 
case however with alternating coefficients. Hence, we found 

\*( y = o)) = \* K ) = j2(- i rn- (9) 

m— 

The ground state in the limit y — > oo is simpler. It is annihilated by Q\ and Q[, which 
act on sites 3j — 2 and 3j — 1 with j — 1 , . . . , n, and thus on the n boxes introduced previously. 
We know that by continuity we must look in the sector with n particles. As the boxes are 
separate, we can try to solve the problem + g? 3 .,_i)|x) = (<4j_2 + °4j-i)lx) = nrs t 

and then take appropriate tensor products. From section [3.2. 1| we already know that this 
problem is solved by a singlet state •■■ 10 — 01 ■■• . Therefore we conclude that the 
y — y oo ground state is given by 

My)) - m = Il(i--2 -4i-i)|o). 

3=1 

4.2.2 Finite coupling 

Let us finally turn to the case of finite coupling. For N — h sites, the ground state 
components are given by 

10100 -y 10010 l + 2y 2 

01010 — (1 + 2y 2 ) 10001 - 2y 2 

01001 l + 2y 2 00101 -y 

Together with the square norm ||*(y)|| 2 = (l+2j/ 2 )(3+8y 2 ). Wehave (<& K \W (y)) = 3(l+2y 2 ) 
and {9s\^[y)) — 3+8y 2 , i.e. both polynomials occurring in the square norm can be obtained 
as projections on asymptotic ground states. Moreover, we see that the state given here is 
even under the parity operation. The fact, that we have definite parity is clear because for 
the present staggering and number of sites the Hamiltonian commutes with parity [H, P] = 0. 
The parity of the ground state at finite y can thus be read off from the asymptotic cases 
discussed in the last section. Upon inspection, we find thus that the ground state for this 
staggering at N = 3n — 1 sites has parity (—1)" (notice that as opposed to the case of 
periodic chains studied in [TH], the singlet state at y — !> oo can be even under parity, namely 
if n is even). 
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Conjecture 4.2 For general n all ground-state components are polynomials with integer 
coefficients and definite parity under y -4 —y. The projection on the kink state ^ is given 
by 



n{ri— 1) 



(* K my)) = (n + l)m(y), m(y)=\^j Qn-iiC 1 ), ( 2 = l + 8y 2 - (10) 

The structure of the wave function is given by 

\*(y))=m(y)\9 K ) + \% (*jf|*)=0, 

and thus considerably simpler than for the kink state of staggering (I), what allows for a 
straightforward application of the iterative algorithm. The square norm factorises into two 
polynomials with positive integer coefficients which coincide with the projections on the kink 
and singlet state respectively. 

\My)\\ 2 = ^my)\*K)(* s my)) 

n + J 

The projection on the singlet state ('5 s\^(y)) yields the unknown polynomials t n (y) introduced 
in section^ 

4.2.3 Densities 

We investigate again the boundary density (y) . As the physics of weak-coupling regime 
is governed by the kink state \^k) we expect the density not to be scale-free at y = 0. 
Indeed, up to n = 8 particles the series expansion is rather compatible with the series 

m . , 1 n-l o n(n- l)(n + 3) 4 3n 4 + 9n 3 + Von 2 + in + 2 fi 

7i+l (n+lj^ (n+l)' 1 (n+1) 4 

4(3n 5 + 12n 4 + 18?i 3 + 20n 2 + lln + 4)y 8 
(n+1) 5 

In the infinite-system limit n -> oo we thus find 
P^(y) = y 4 -3y 6 + 12y 8 + ... 

what matches the first terms of the series expansion at y = of the corresponding results for 
N = 3n sites given in Q. Based on physical grounds, we conjecture that our the boundary 
density considered here converges to the expressions given there. 

Conversely, the boundary density in the y > 1 regime displays at scale-free behaviour 
around the point y -4 oo. From the stable part of the series expansion we conjecture the 
expression 

(i), , = 1 - V + 8;/ 4 - 8y(y 2 - l) 3 / 2 
Po ° [V> 16y 2 



5 Conclusion 

In conclusion, we presented several conjectures on the ground-state properties of a staggered 
one-dimensional supersymmetric fermion model with open boundary conditions, based on 
the exact diagonalisation of small systems. In particular, we conjecture a relation between 
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certain ground state components of the model, and a class of polynomials which appeared 
in the literature on the XYZ spin chain along a special line of couplings. The very same 
polynomials occur in the square norm of the corresponding ground state vectors. Moreover, 
exploiting the scale-free property we determined exact expressions for local site-occupation 
probabilities for infinite systems. 

Given these observations it is quite natural to ask if there are open boundary conditions 
for the XYZ spin chain, possibly amongst the integrable ones [27], which still lead to exact 
zero-energy ground states with similar properties to the ones studied here. This is the case 
for the XXZ chain as was shown in [S] , yet it might be more difficult off criticality because 
one can show that the underlying supersymmetry in the spin chains does not survive this 
deformation |25j . Another interesting open problem concerns the choice of the coupling 
constants Aj-: the results presented here relate to a staggered pattern. In fact, it turns out 
that some of the conjectured sum rules survive for generic choices. In fact, one may show 
that these follow from the supersymmetry alone. This problem will be addressed in [2"B] . 
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A Polynomials 

In this appendix, we list various polynomials encountered in the text pk{w) and qk{w). They 
can also be found in [20] • The first explicit cases of the pk(w) polynomials are 

p-i(w) = 1 - Aw + 22w 2 - 70w 3 + 245w 4 - 680w 5 + 1940w 6 - 3260w 7 + 5915w 8 



p 2 (w) = l + 2w + 7w 2 + 10w 3 + 21w 4 + 12w 5 + llw 6 , 

ps(w) = 1 + 3w + I5w 2 + 35u> 3 + 105w 4 + 195w 5 + 435w 6 + 555w 7 + 840w 8 
+710w 9 + 738w 10 + 294W 11 + 170w 12 , 

Pi{w) = 1 + Aw + 26w 2 + 82w 3 + 319w 4 + 840w 5 + 2488w 6 + 5572w 7 + 13524w 8 
+24920w 9 + 48776w 10 + 72800U; 11 + 114716w 12 + 135464w 13 + 169536w 
+148972w; 15 + 141835w 16 + 85044w 17 + 58406w 18 + 17822w 19 + 7429w 2C 



p-i(w) 
px(w) 



-6500w 9 + 7878w 10 - 3822W 11 + 2431w 12 , 
1 - 3w + 12w 2 - 30u> 3 + 81w 4 - 63w 5 + 66w 6 , 
1 - 2w + 5w 2 , 

1, 
1, 

1 + w + 2w 2 , 
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The first explicit cases of the qk(w) polynomials are 

q-4(w) = 1 + Sw 2 + 35w 4 + 560w 6 + 1215w 8 + 1848w; 10 + 429w 12 , 

q- 3 (w) = 1 + 3w 2 + 39u; 4 + 2lw 6 , 

q- 2 {w) = l + 3w 2 , 

q-i(w) = 1, 

<7oO) = 1, 

qi(w) = l + 3w 2 , 

q 2 (w) = 1 + 8w 2 + 29w 4 + 26w 6 , 

q 3 {w) = l + 15w 2 + 112w 4 + 518w; 6 + 1257w 8 + 1547w 10 + 646w 12 , 
q 4 {w) = 1 + 24w 2 + 291w 4 + 2338w 6 + 13524w 8 + 54474w 10 + 150472w 12 + 276678w 
+312195w 16 + 192694w 18 + 45885w 20 . 
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